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Abstract: The single inclusive spectrum of charged particles with transverse momenta 
Pj, = 3 - 150 GeV/c measured at midrapidity by the CDF experiment in proton-antiproton 
(p-p) collisions at ^/s = 1.96 TeV is compared to next-to-leading order (NLO) perturbative 
QCD calculations using the most recent parametrizations of the parton distributions and 
parton-to-hadron fragmentation functions. Above Pj, ~ 20 GeV/c, there is a very sizeable 
disagreement of the Tevatron data compared to the NLO predictions and to x^-scaling ex- 
pectations, suggesting a problem in the experimental data. We also present the predictions 
for the Pj,-differential charged hadron spectra and the associated theoretical uncertainties 
for proton-proton (p-p) collisions at LHC energies {^/s = 0.9 - 14 TeV). Two procedures 
to estimate the charged hadron spectra at LHC heavy-ion collision energies (y^ = 2.76, 
5.5 TeV) from p-p measurements are suggested. 
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1. Introduction 

Hadron production at large transverse momenta {pj. ^ Aqcd ~ 0.2 GeV) in hadronic 
interactions originates from the fragmentation of the hard scattered partons produced in 
the collision. The presence of a hard scale in the process allows one to employ the powerful 
theoretical machinery of collinear factorisation [1] to compute the corresponding production 
cross sections. High-p^ hadron cross sections can be thus obtained as a convolution of (i) 
long-distance universal pieces representing the structure of the initial hadrons (parton 
distribution functions, PDFs) as well as the fragmentation of a final-state quark or gluon 
into the observed hadron (fragmentation functions, FFs), and (ii) short-distance parts that 
describe the hard partonic interactions calculable as a perturbative expansion in terms 
of the strong running coupling . The measurement of high-p^ hadroproduction in p-p 
and p-p collisions provides, thus, a valuable testing ground of the perturbative regime 
of Quantum Chromodynamics (pQCD) and of the non-perturbative objects (PDFs, FFs) 
needed to compute a large variety of cross sections at hadronic colliders. 

Theoretically, lowest-order (LO) calculations of the inclusive hadron cross sections were 
performed in the late 70s [2], later improved at next-to-leading order (NLO) [3-5] and 
more recently at next-to-leading-log (including soft gluon resummation) [6, 7] accuracies. 
The latest phenomenological developments in this field have focused on constraints of the 
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proton PDFs (in particular the polarised ones [7, 8] ) , on improvements of the parton-to- 
hadron (in particular, gluon-to-hadron) FFs [9], as well as on baseline measurements of 
relevance for high-energy heavy-ions collisions [10]. On the experimental side, inclusive 
unidentified charged hadron production - i.e. pp,pp — t- h^X, where = {h^ + h^) is 
effectively the sum of pions (about 60% of all hadrons), kaons (about 20% of the total) and 
protons (about 10% of all hadrons) and their antiparticles - have been measured above 
^ 1 GeV/c at the ISR {^/s = 31, 44, 63 GeV) [11], at RHIC {^/s = 200 GeV) [12], SppS 
i^/s = 0.2, 0.5, 0.9 TeV) [13], and Tevatron (^i = 0.63, 1.8, 1.96 TeV) [14-16] energies. 
Except at Tevatron, the rest of measurements are unfortunately in a moderate Pj, range 
{Pj, ~ 12 GeV/c at most). The latest comparisons of the available charged hadron spectra, 
at RHIC energies [8,12], to NLO calculations [6] show a good data-theory agreement above 
Pj, 1.5 GeV/c for central and forward rapidities [10]. 

In this paper, we compare NLO pQCD calculations to the latest charged particle spec- 
trum measured at Tevatron and we present predictions with their theoretical uncertainties 
for the high-pj, hadron spectra expected at LHC energies. The motivation is two-fold. 
First, the most recent CDF charged particle spectrum [16] covers a very large Pj, range, 
up to Pj. = 150 GeV/c where pQCD predictions are reliable and can be confronted to 
the data. Similarly, comparable "minimum bias" measurements are expected to be avail- 
able in the early running of the LHC [17]. CMS has already measured a first, yet mostly 
low-pj,, charged hadron spectrum at ^/s = 2.36 TeV [18]. Secondly, at the LHC, a p-p 
reference hadron spectrum will be needed at the same centre-of-mass (cm.) energy as 
that of heavy-ion (Pb-Pb) collisions to study the high-p^ suppression observed in nucleus- 
nucleus reactions at RHIC [19]. Since the Pb-Pb results will be nominally obtained at 
^/s^ = 5.5 TeV a pQCD-based interpolation between the results recorded at Tevatron 
{^/s = 1.96 TeV) and during the first LHC p-p run (y^ = 7 TeV) will be needed. 

The paper is organized as follows. In Section |2| we succinctly remind the theoretical 
framework of our study based on the next-to-leading-order pQCD Monte Carlo (MC) code 
INCNLO [20]. In Section |^ we compare the charged particle spectra measured at mid- 
rapidity in p-p collisions at y/s = 1.96 TeV [16] to the NLO calculations incnlo and to the 
LO parton shower MC pythia, as well as to simple x^-scaling expectations. We find that 
the maximum theoretical uncertainties of the NLO prediction - associated to the PDF, FF 
and scale variations added in quadrature - are ±40%. For increasing transverse momenta, 
the data is a factor up to 3 orders of magnitude larger than the perturbative predictions. 
We conclude that above Pj. « 20 GeV/c, there is no possibility to accommodate the data- 
theory discrepancy even accounting for possible additional contributions to the charged 
particle yield coming e.g. from heavy-quarks or vector-boson (plus jet) production. The 
fact that the parent jet -differential spectrum is, on the contrary, well reproduced by NLO 
calculations and that the single particle data also violate simple x^-scaling expectations, 
suggest a problem in the experimental results at the highest Pj, values. Finally, in Section ^ 
we present the charged hadron spectra and associated uncertainties predicted by incnlo in 
p-p collisions in the range of energies covered by the LHC {^/s = 0.9 - 14 TeV), and propose 
two methods to determine the p-p spectra at intermediate LHC energies of relevance for 
heavy-ion running. 
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2. Hadroproduction in factorised pQCD 



The inclusive cross section for the production of a single hadron, differential in transverse 
momentum p^, and rapidity y, takes the following form at next-to-leading order in [21]: 



/dz 
dxi dx2 Fi/p{xi,Hp) Fj/p{x2,Hp) ^ D^{z,fi^ 



da _ _ 



+ ( TTl ) ^ij,k{^^R■:^^F■'^^FF) 



2tt J dp^dy \ 27r 



(2.1) 



Here i^j/p(xi, /i^) are the PDFs of the incoming protons p at parton momentum fraction 
D^{z, Upp) are the parton-to-hadron FFs describing the transition of the parton k into 
an unidentified hadron h carrying a fraction z of its momentum, daij^k/dp^dy is the Born 
cross section of the subprocess i -|- j — )• /c -|- X, and Kij^k is the corresponding higher-order 
term (the full kinematic dependence is omitted for clarity). In this paper, we use the 
INCNLO programme [20] to compute the cross sections, supplemented with various PDFs 
and FFs sets (see below). The truncation of the perturbative series at next-to- leading order 
accuracy in introduces an artificial dependence with magnitude O (og^), of the cross 
section on initial-state (np) and final-state {fipp) factorization scales, as well as on the 
renormalization scale . The choice of scales is arbitrary but the standard procedure is to 
choose a value around the natural physical scale of the hard scattering process, here given 
by the p^, of the produced hadron. We consider below scale variations fi^^, n^, fipp = npp, 
with K = 0.5 - 2 to gauge the theoretical uncertainty linked to the neglected higher-order 
terms. Hereafter, whenever the scales /i^, ^p and ^pp, are given a common value, the 
latter is denoted jj,. 



The two non-perturbative inputs of Eq. (2T) are the parton densities and the fragmen- 



tation functions. The former are mostly obtained from global-fit analyses of HERA proton 
structure function data, the latter from hadron production results in e^e^ collisions. We 
use here the three latest PDFs parametrisations available: CTEQ6.6 [22], MSTW08 [23] 
and NNPDF1.2 [24], included in the lhapdf (version 5.7.1) package [25], which take into 
account the most up-to-date data from deep-inelastic lepton-proton scattering and hadronic 
collisions as well as various theoretical improvements. For the fragmentation functions into 
hadrons, we use and compare the three more recent FF sets available: DSS [26], AKK08 [27] 
and HKNS [28], which, except for the latter, include for the first time also hadron-hadron 
collision data in their global analyses. These new FF fits cover a larger z range and are 
more sensitive to the gluon fragmentation which dominates high-p^ hadron production in 
p-p collisions [9]. 

For transverse momenta close to the phase space boundary where the pT of the hadron 
is about half of the partonic centre-of-mass energy {x^ = 2p^/ y/s ^ 0.1 — 1), the coefficients 
of the perturbative expansion are enhanced by extra powers of logarithmic terms of the 
form ag"'ln^"~^(l — Xp)/{1 — Xp) [29]. Resummation to all orders of such "threshold" 
terms - which appear when the initial partons have just enough energy to produce the 
high-pj, hadron - have been carried out at next-to-leading logarithmic (NLL) accuracy 
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in [6,7]. Interestingly, the NLL results provide a much reduced scale-dependence than 
the NLO approximation. The presently used fixed-order calculations (incnlo) do not 
include threshold resummations but their effect in the final spectrum is expected to be less 
important since the typical charged hadron Pj, range covered by the Tevatron and LHC 
experiments, = 2py/y^ ~ 10~^ — 10""*^, is far away from the region where such effects 
start to play a role. 



3. Tevatron data versus perturbative QCD 

In this Section we compare the high-p^, charged particle spectrum measured by the CDF 
collaboration [16] in the pseudorapidity range |r/| < 1 to the predictions of incnlo [20] 
and PYTHIA [30] MCs and to simple perturbative expectations based on Xj,-scaling [31]. 
The measured spectrum covers the range pT = 0.4 - 150 GeV/c, but a comparison to 
pQCD predictions is only meaningful at high enough pT] therefore we impose a minimal 
cut of = 3 GeV/c. For the NLO analysis, we study separately the effects on the 
spectrum of varying in the calculations the three theoretical scales (/i = Pj,/2,pr^,2pj,), 
PDFs (MSTW08, CTEQ6.6 and NNPDF1.2) and FFs (AKK08, DSS and HKNS). We use 
PYTHIA to determine possible extra contributions to the measured high-p^, tracks spectrum 
coming from heavy-quark fragmentation as well as from real and virtual vector-boson 
production, either single-inclusive or in association with a jet. 

3.1 Data versus incnlo 

The measured CDF charged particle p-p single inclusive distribution is compared to the 
INCNLO predictions for charged hadrons in Fig. ||. First, we note that the measured pri- 
mary track spectrum is not corrected for contributions from charged particles other than 
hadrons. Possible contamination from stable leptons (electrons, positrons and muons) 
are not in principle removed from the measured spectrum. As we discuss a posteriori in 



Section 3.2, those amount however only to a small fraction (a few percents) of the total 
charged particle tracks coming from quark and gluon jet fragmentation according to our 
PYTHIA simulations. The incnlo prediction shown in Fig. || is that which best fits the 
(low Pj, range of) the experimental results. We see that below p^, w 20 GeV/c data and 
theory agree well for the choice of scales fx = 2p^, CTEQ6.6 parton densities, and AKK08 
parton-to-hadron fragmentation functions. Above this p^ value, the CDF spectrum starts 
to rapidly deviate from the theoretical predictions. At the highest transverse momenta the 
data is up to a factor 800 above the NLO calculations. A very conservative quadratic sum 
of all uncertainties discussed hereafter - amounting to ±30% for the scales, ±10% for the 
PDFs, and ±25% for the FFs choices - would result in a maximum theoretical uncertainty 
of ±40% in the yields (dashed lines around the data/theory ratio). 

Scale uncertainty: First, in the incnlo calculations we have fixed the PDFs and FFs 
to the CTEQ6.6 and DSS sets respectively^, and computed the corresponding spectra 



^The choice is in principle arbitrary, other PDF and FF combinations yield similar results for the scale 
dependence. 
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Figure 1: Top: Comparison of the charged particle spectrum measured by CDF in p-p cohisions 
at = 1.96 TeV [16] to NLO pQCD predictions with PDFs fixed to CTEQ6.6, scales to = 2p^, 
and FFs to AKK08. Bottom: Corresponding ratio of CDF data over theory. The dashed lines 
indicate the maximum ±40% theoretical uncertainty of the calculations (see text). 
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INCNLO: PDF=CTEQ6.6, FF=DSS 
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pp^h- X,\/s= 1.96 TeV, |ri|<1 



4 5 6 7 8 10 



20 30 40 50 



100 200 
(GeV/c) 



C 0.8- 

"lo.e: 

0) 

O 0.4 

c 

D 0.2- 



INCNLO:n=p^, FF=DSS 

□ CTEQ6.6 vs. MSTW08 

□ CTEQ6.6 vs. NNPDF1.2 

□ MSTW08 vs. NNPDF1.2 



Q 
Q. 




0.2 
-0.4 
-0.6 
-0.8 
-1 

E 0.8 

2 0.6 
(D 

O 0.4 
C 

3 0.2 

t 

-0.2 



pp^ii-X,xJs= 1.96 TeV, |ri|<1 



3 4 5 6 7 8 10 



20 30 40 50 100 200 

(GeV/c) 



-1 





INCNLO: ^= p^, PDF=CTEQ6.6 






□ HKN vs. AKK08 






□ DSS vs. AKK08 






□ DSS vs. HKN 










pp^h-X,Njs= 1.96 TeV, |ri|<1 











3 4 5 6 78 10 



20 30 40 50 100 200 

p^ (GeV/c) 



Figure 2: Fractional differences between the incnlo charged hadron spectra in p-p at 
y/s — 1.96 TeV for varying scales fii, PDF and FF. Top: Scale uncertainty obtained for fixed 
PDF (CTEQ6.6) and FF (DSS) varying all three scales within = - (the dashed lines 

indicate a ±30% uncertainty). Middle: PDF uncertainty obtained for fixed /i = and FF (DSS) 
with three PDFs: CTEQ6.6, MSTW08, NNPDF1.2 (the dashed hues indicate ±10% differences). 
Bottom: FF uncertainty obtained for fixed scales {fi — Pj.) and PDF (CTEQ6.6) with three FFs: 
AKK08, DSS, HKNS (the dashed lines represent&t25%). 



setting the three theoretical scales to three different values iJipj^pp = Pj,/2,p,^,2p^ in all 
possible 27 combinations. The corresponding range of predictions is shown in Fig. ^ (top) 
where we plot a shaded band covering the whole range of fractional differences between the 
spectra obtained for any choice of scales. The "closest-to-the-average" spectrum is obtained 
setting all scales to ^ = p^.. The largest (resp. lowest) charged hadron yield predictions are 
obtained with mostly all three scale values set to /Uj = Pj./2 (resp. = 2p^). At low 
the scale uncertainty is quite large (indicating as expected larger higher-order corrections) 
but otherwise above = 10 GeV/c it stays roughly constant at around ±30% up to the 
highest momenta considered (dashed lines in the figure). 

PDF uncertainty: Second, in the middle panel of Fig. ^ we show as a function oip^, the 
theoretical uncertainty associated to the PDF choice. It has been obtained with INCNLO 
comparing the fractional differences between the single charged hadron spectra at ^/s = 1.96 
TeV for fixed scales (fJ- = Pj,) and FF (DSS) and three different PDFs. The dashed bands 
plotted cover the range of maximum relative differences in the theoretical spectra obtained 
with MSTW08, CTEQ6.6 and NNPDF1.2. Those differences are small, below 10% and 
mostly py-independent (dashed lines). 

FF uncertainty: Last, we have used incnlo complemented with the three latest FFs 
available in the market: AKK08, DSS and HKNS, to compute the charged hadron spectrum 
for the CDF kinematics, with scales (/i = Pj,) and PDFs (CTEQ6.6) fixed. The main 
differences between FF sets concern the fractional vr^, and p/p compositions as a 
function of Pj,. Yet, the total^ hadron yield predicted by the three FFs for p-p at 1.96 TeV 
is quite similar as can be seen in the bottom panel of Fig. ^ where we plot the relative 
differences between the spectra computed for varying FFs. The maximum theoretical 
uncertainty linked to the FF choice amounts to about ±25% of the charged hadron yield 
at any p^ (dashed lines), although the differences between FFs appear to be smaller at 
intermediate hadron p^, ~ 10 - 25 GeV/c. 

A conservative quadratic sum of the fractional uncertainties linked to the NLO scales, 
PDFs and FFs choices results in a total theoretical uncertainty of order ±40% whereas 
the maximum difference between the data and the calculations amounts to much larger 
factors, up to 0(10^) at p^ > 100 GeV/c (Fig. |lD. 

3.2 Data versus pythia 

Given the large discrepancy between the experimental and NLO predictions for the charged 
particle spectrum at high-p^ one may wonder whether other charged particles - apart 
from TT^,K^,p/p coming from the fragmentation of quarks and gluons - may contribute 
in any way to the experimentally measured distribution beyond ~ 20 GeV/c. A first 
possibility that we have considered is whether other charged products from charm and 
bottom jets (with relative increasing importance at large transverse momentum) play any 

^As a cross check, for all FF sets we have confirmed that the NLO spectrum obtained from the sum of 
the spectra individually obtained with the pions, kaons and protons FFs is indeed equal to the one obtained 
with the non-identified charged hadron FFs. 
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role. Although the inclusive hadron FFs used in our NLO calculations contain all pion, 
kaons and (anti)protons issuing from light- as well as heavy-quark fragmentation, charm 
and bottom hadrons decay also into charged leptons^ which are not included in INCNLO. 
Thus as a independent theoretical check, we have computed the inclusive yield of all charged 
particles with the pythia MC (v6.420) [30] with the D6T [33] tuning"^ in the "minimum 
bias" and QCD-jets modes (msel = 1 with \ow-pj, production, ISUB = 95, switched on to 
correctly simulate the low-p^ region). The chosen processes produce not only light-quarks 
and gluons but also heavy-quarks^ including flavour excitation, Qg — t- Qg, and gluon 
splitting, g — )• QQ. 

To obtain enough statistics at high-p^ , we have run with up to 12 different ranges for 
the minimum and maximum parton momenta in the 2— >-2 scatterings {pj, = - 10, 10 - 15, 
15 - 20, 20 - 50 GeV/c,. . .and Pj, > 470 GeV/c) weighted by their corresponding cross 
sections. We have then explicitly separated the contributions coming from the fragmenta- 
tion of high-py light-flavours {u, d, s and gluon) from those coming from charm and bottom 
quarks. As done in CDF, we take all charged particles^ exactly as defined in their analysis 
(i.e. all primary particles with mean lifetimes r > 0.3 10^^*^ s and the decay products 
of those particles with shorter r). The results of our studies are shown in Fig. ^. The 
inclusive charged products of c-quark and 5-quark fragmentation represent a very small 
(less than 5%) fraction of the total yield of particles measured at high-p^, by CDF. 

A second possibility that we have explored is whether the charged products of real 
and virtual vector-boson (7^*\ W^, and Z^) production - either single-inclusive or in as- 
sociation with a jet - which start to play a role at increasing transverse momenta, could 
partially account for the data-theory discrepancy. We have run prompt photon produc- 
tion in PYTHIA including the Born-level 7-jet Compton and annihilation diagrams (ISUB = 
29,14 respectively). The W^, Z° and DY production (iSUB = 1,2,15,16,30,31), includes 
single-inclusive (2 — t- 1) as well as double-inclusive (2 — t- 2) VF^-,Z'^-,DY-jet channels with 
Pt > 20 GeV/c. The and contributions produce a (local) Jacobian peak in the 
charged-particle p^ distribution at about half the vector-boson mass, Pj, ^ 40 GeV/c. All 
those contributions, shown added up in Fig. ^, increase the charged particle yield by up 
to 10% in the range above Pj, w 40 GeV/c. This number is consistent with a simple order 
of magnitude estimate based on the ratio of electroweak and strong coupling constants 
valid when p^ > Mw/2: {a^^/usf = (0.034/0.12)^ = O (lO-^). Clearly, those processes 
contribute little to the total yield of charged particles and therefore cannot justify the 
observed large discrepancy between data and theory. 



^As a cross check, we have confirmed that the pythia spectrum of single leptons from c and h production 
agrees relatively well (within a factor of two) with more involved fixed-order NLL calculations [32]. 

*Tune D6T uses the CTEQ6LL PDF and describes the underlying event and the Drell-Yan data at 
Tevatron. 

^PYTHIA settings: PARP(91)=2.1 GeV/c (intrinsic kr), PMAS(4, 1)=1 . 5 GeV/c^ (m^ mass), 
PMAS(5,1)=4.8 GeV/c^ (ms mass), MSTP(33)=1 (Ti'-factor). Alternative running of standalone heavy- 
quark production (with MSEL = 4 and 5) would require A'-factors of 2 - 4 in order to reproduce the heavy 
flavour spectra measured at various colliders [34] . 

^PYTHIA settings: MSTJ(22)=2, PARJ(71)=10. 
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Figure 3: Top: Comparison of the CDF data (symbols) to PYTHIA (v6.420, D6T tuning) 
distribution of charged particles in p-p collisions at y/s = 1.96 TeV coming from the fragmentation 
of (i) light-quarks and gluons and (ii) heavy-quarks, and (iii) from processes involving the production 
of vector bosons (7^*^ , W^, Z^). Bottom: Ratio of the CDF data to the sum of all pythia charged 
particle contributions. 
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3.3 Data versus a^^-scaling 



A robust pQCD prediction for hard processes A B ^ C X in hadronic collisions is the 
power-law scaling of the inclusive invariant cross section, 

E Sa/d^p = F(x,)/p?"-^^ = . (3.1) 

In the original parton model the power-law fall-off of the spectrum is simply n = 4 since 
the underlying 2—7-2 subprocess amplitude for point-like partons is scale invariant. In 
QCD, small scaling violations appear due to the running of and the evolution of PDFs 
and FFs. At midrapidity and at fixed = 10 GeV/c, the power-law exponent computed 
at NLO accuracy increases slowly from n ~ 5 at small values of Xj, {x^ = 10"'^) up to 
n ~ 6 at Xj, = 0.5, with a very small dependence on the specific hadron species [35]. 




Figure 4: Compiled charged particle cross-sections measured in p-p collisions at five different cm. 
energies from 0.2 TeV to 1.96 TeV plotted as a function of (left) and as a function of (right) 
scaled with an effective common exponent of rt = 5.5 (see text). 



Except for the latest CDF data, the theoretical expectation, Eq. (3J.), is indeed well 
fulfilled by the experimental charged particle spectra measured so far in p-p collisions'' at 
different centre-of-mass energies at the CERN SppS {y/s = 0.2, 0.5, 0.9 TeV) [13] and 
Tevatron (y^ = 0.63, 1.8, 1.96 TeV) [14,16] colliders. All the p^, spectra feature power-law 
behaviours above Pj, ~ 2 GeV/c (the higher the cm. energy the smaller the exponent 
of the fall-off, see Fig. Q left). Following the expectation Eq. (p.l|), in order to extract a 
common n value from these different data sets, the measurements are plotted as a function 



factor of 1/2 is applied hereafter to tlie CDF Run II spectrum as they measured {h^ + h ) instead 
of the average 0.5 x + h~) for all other measurements. 
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of Xj,, multiplied by ^/s"' and fitted with the following 3-parameter functional form 



l-n 



y=0 



PO • [1 + ixT/pi)r ■ 



(3.2) 



In the data fitting, a minimum px of 2 GeV/c is applied to exclude the region where soft 
particle production (which does not follow -scaling) is dominant, which is consistent 
with what is used in [35]. We also exclude the CDF Run-II data from the global fit since, 
as we see a posteriori, there is no possibility to get an agreement with the lower energy 
measurements. With the obtained {pj}-parameters, using n^^o ~ 5 as a guidance, the 
exponent n is varied from 4 to 7 in incremental steps in order to minimize the following 
function with MINUIT [36] 



ndat 



X in,{Pi}) 



y=o 



poii+(£TVpi)Fi H 2 



(3.3) 



where aj are the quadratic sum of the statistical and systematic experimental uncertainties. 
In Fig. ^ (right) we show the experimental charged particle spectra scaled by y/s^ as a 
function of Xj, with the best value of n obtained from the fit, n = 5.5. We note that all 
measurements spanning a range of one order-of-magnitude in centre-of-mass energies follow 
a universal curve after rescaling up to the highest x^, ~ 0.03 measured at lower energy. A 
deviation of the CDF measurement at ^/s = 1.96 TeV from the trend established by the 
lower measurements is prominent above Xj. ~ 0.03 {pj. ~ 30 GeV/c). 

In order to better assess the amount of deviation of the CDF data to the Xj, -scaling 
expectation we show in Fig. |5| (top) the xy-scaled fit obtained from all lower-energy data 
extrapolated to an expected Pj, spectrum at ^/s = 1.96 TeV compared with the CDF Run II 
measurement and with the PYTHIA prediction shown in Fig. ^ As observed before, beyond 
Pj, ~ 20 GeV/c the latest CDF data clearly fail to follow the Xj,-scaling expectation fulfilled 
by the rest of charged hadron measurements. 



3.4 Discussion 

The fact that no combination of PDF, FF and/or theoretical scales in the NLO calculations 
is able to reproduce the Tevatron experimental data above p^, ^ 20 GeV/c by such a large 
factor is totally unexpected. Indeed, similar calculations based on fastnlo/nlojet-1-+ [37, 
38], reproduce perfectly well the inclusive jet spectrum measured in p-p collisions at 
^/s = 1.96 TeV in the range p^^* 50 - 600 GeV/c. The shape and magnitude of the 
CDF jet measurement is well reproduced using CTEQ6.1M PDFs and renormalization and 
factorization scales^ set to /U^ = //^ = p!^^/2 [39]. Likewise, the D0 jet measurement agrees 
well with the same NLO predictions with CTEQ6.5M parton densities and /i^ = /U^ = p!^^ 
scales [40]. Variations of PDF and/or scales in the jet calculations, result in differences 
typically of order 10% - 15% for both measurements [39,40]. 

*Note that in the inclusive jet calculation there is one scale less, the fragmentation one fip. 
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Figure 5: Top: Inclusive CDF Run-II charged liadron spectrum (filled circles) compared to the 
a;T-scaled extrapolation of lower energy p-p data (blue solid line), as well as to the total pythia 
prediction (dashed red line) of Fig. ^. Bottom: Corresponding ratio of CDF data over x^-scaling 
expectation. 
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Given the agreement between the jet data and the fixed-order calculations^, it is some- 
how difficult to conceive a strong disagreement in the hadron production channel since the 
single high-pj, charged particle spectrum is dominated by leading hadrons carrying out a 
large fraction, {z) ~ 0.6-0.7, of the parent parton energy. Quite naively, the distribution of 
charged hadrons above ~ {z) p-!^^ > 30 GeV/c should be also perfectly consistent with 
the theoretical predictions within the additional uncertainty introduced by the fragmenta- 
tion functions which is at most of the order ±25% as seen in Fig. I (bottom). Let us give 
a rough estimate of the expected invariant hadron production cross section based on the 
jet data. At leading order accuracy and assuming that one partonic channel dominates jet 
and hadron production (e.g. uu — )• uu scattering at large Xj,), the hadron production cross 
section Eq. is roughly given by 



d3p ~ ^ If J <^-''-' ^ ^ ^ d3p 



with z ~ 0.7, the typical range Az ~ 10 ^ which contributes to the hadron production 
cross section, and D^~^^ {z,Pj,) ~ 10^^ the u-quark-to-hadron FF (see e.g. [42]). Using 
the D0 jet measurement of d'^a^^ydp^dy ~ 200 pb/(GeV/c) at p^ = 100/z ~ 140 GeV/c 
and \y\ < 1 [40], one gets for the hadron production cross section a value of d^a'^/d^p ~ 
2 X 10~^^ mb/(GeV/c)^ which is the right order of magnitude estimate expected in QCD 
(see the LO pythia curve in Fig. The inconsistency between jet and the CDF large-p^ 
hadron spectrum is also discussed in detail in [43,44]. 

Of course the above argument relies on the factorization assumption that large-p^ 
hadron production can be expressed as a convolution of hard matrix elements with parton- 
to-hadron fragmentation functions. Should the FFs - mostly based on fits of e^e~ data - 
be non-universal one could imagine that the recent CDF measurement actually reflects 
dramatic modiflcations of fragmentation functions in hadronic collisions. This however 
seems unlikely given the success of the DSS [26] and AKK08 [27] global fits of fragmenta- 
tion functions which consistently use e~^e~ data together with RHIC measurements in p-p 
collisions at y/s = 200 GeV. It is in particular unclear how possible factorization breaking 
effects could enhance hadron production cross sections by up to 3 orders of magnitude at 
transverse momenta as large as p^ = 150 GeV/c. 

There exist even more general arguments why the large-p^, CDF data cannot be under- 
stood as coming from hadron production in perturbative QCD. As discussed in the previous 



section (Sect. |3.3D , the recent CDF spectrum departs from the x^-scaling behaviour ob- 
served in the lower-energy data which can all be described assuming a scaling exponent 
n ~ 5.5. We find on the contrary that the scaling exponent obtained from the comparison 
of the large-py CDF measurement with the UAl data at ^/s = 200 GeV is roughly^*^ n ~ 4- 
4.7. This value is extremely close to what is expected in the conformal limit (n = 4), i.e. 
assuming no scaling violations at all in QCD. It is in particular smaller than the exponents 



^Prompt photon data at the Tevatron are also very well reproduced at large by NLO pQCD calcu- 
lations [41]. 

^"The precise value is difficult to obtain since the a;^ -spectra at the two cm. energies have a different 
slope, already indicating a non-conventional behaviour in one of the two data sets. 
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expected for jet and prompt photon production [35], despite the fact that scahng violations 
are expected to be stronger in the hadron production channel because of the additional 
fragmentation process. What is more, the scaling exponent n obtained at fixed Xj, reflects 
the -dependence of the hard partonic cross section a ~ Because of the fast variation 
of the parton densities with x^,, the p^-slope, a, of the invariant production cross section 
at fixed is expected to be somehow larger than the scale dependence of the partonic 
cross section, i.e. a > n. Surprisingly the value a obtained from a fit of CDF data alone 
above Pj, = 22 GeV/c is as small as a ~ 3.9, that is smaller than the combined UAl-CDF 
scaling exponent n ~ 4-4.7 (and even lower than the smallest scaling exponent n = 4!). 
This clearly indicates that it is not possible to describe the Run-II CDF measurement as 
coming from hadron production in perturbative QCD. Hence factorization breaking effects 
in the fragmentation channel cannot be at the origin of the present discrepancy between 
data and NLO theory. 

We conclude from this Section that the facts that (i) the NLO calculations largely 
fail to reproduce the measured single-hadron spectrum at large p^, while reproducing well 
the single jet p^, -differential cross sections, and (ii) that the measurement violates simple 
phenomenological expectations such as -scaling confirmed empirically in all hadronic col- 
lisions so far, point to a possible experimental problem in the data above Pj, ~ 20 GeV/c 
- or from unknown sources of charged particles - rather than from a sudden breakdown of 
QCD perturbation theory in the hadron production channel. 



After we finished this work, other analyses appeared [43-46] that point out to the same 
discrepancy between the CDF data [16] and NLO calculations. In [45], Albino, Kniehl and 
Kramer point out the disagreement between data and NLO theory and question the validity 
of factorization theorems for large-pj, hadron production, a possibility which we exclude 
(see discussion above). In [43] it is shown on general grounds that the spectrum measured 
by CDF is inconsistent with existing Tevatron data on the inclusive jet production cross 
section and the distribution of hadrons inside jets (a similar argument is given in [44]). 
This observation allows the authors to exclude, as well, the breakdown of factorization as a 
possible explanation of the data. They also conclude that new physics scenarios explaining 
the CDF excess are unlikely, yet they cannot be fully eliminated. Finally, it has been 
claimed in [46] that weak boson decays into hadrons might explain the CDF data. This 
possibility is however excluded as shown in Section 3.2, either from the detailed pythia 
calculations or from the order of magnitude estimate^ ^. 



4. Inclusive charged hadron spectra at the LHC 

In this last Section of the paper we present first the incnlo predictions for the charged 

^^We believe that the calculations in [46] is incorrect partly because of the use of fragmentation functions 
which are two orders of magnitude larger than the usual fits from e"'"e~ data, at large z. We also note that 
the Jacobian peak in the p^-spectrum is located at ~ mw instead of ~ mw/^. Finally it is difficult to 
conceive why the invariant cross section scales as ~ s/p^ instead of conventional behaviour ~ (we 
thank S. Brodsky for pointing this out), which might lead to another two order of magnitude, O (s/p^), 
overestimate in [46]. 
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hadron p^, -differential cross sections at mid-rapidity (|?7| < 1) in p-p collisions in the range of 
CERN LHC energies {y/s = 0.9 - 14 TeV) including their expected theoretical uncertainties. 
Second, we discuss two interpolation methods, based on pQCD-ratios and Xj,-scaling, that 
can be used to obtain a baseline charged hadron p-p spectrum at intermediate LHC energies 
(^/s = 2.76, 5.5 TeV) needed to compare against similar measurements to be carried out 
in Pb-Pb collisions. 

4.1 INCNLO predictions 

Figure ^ (left) shows the incnlo p^, -differential cross sections in pp — )• h^X at six differ- 
ent cm. energies (expected to be) reached at the LHC at various stages of the collider 
programme. The spectra have been obtained with CTEQ6.6 PDFs, DSS FFs and the- 
oretical scales set to fi = Pj,. Whereas at Pj, below about 10 GeV/c all calculations 
converge, with increasing cm. energies (and thus phase-space for hard parton-parton scat- 
terings), the spectra become increasingly flatter. For example, the charged hadron yield 
at Pj, w 100 GeV/c, increases by a factor of 10 between ^/s = 2.76 TeV and ^/s = 7 TeV 
and yet by another factor of 5 between the latter and the top LHC energy of 14 TeV. As 
shown in Fig. ^ (right), for this particular choice of scales/PDFs/FFs, a common power-law 
exponent of n = 4.9 allows one to scale all NLO spectra in the range ^/s = 0.9 - 14 TeV 
to a universal curve, using the Xj, prescription given by Eq. ( |3.lD . This value is slightly 



smaller than what has been obtained in Section 3.3 from Tevatron and SppS measurements 



(n = 5.5), indicating as expected smaller scaling violations at larger cm. energy [35]. 
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Figure 6: Charged hadron spectra in p-p collisions at ^/s = 0.9, 2.76, 5.5, 7, 10 and 14 TeV 
predicted by NLO pQCD calculations with CTEQ6.6 parton distribution functions, DSS fragmen- 
tation functions, and scales set to fi — p^^: -differential (left) and z^-scaled with exponent n = 4.9 
(right). 
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Figure 7: Fractional differences between the incnlo charged hadron spectra in p-p at y/s = 7 TeV 
for varying scales Hi, PDF and FF. Top: Scale uncertainty obtained for fixed PDF (CTEQ6.6) 
and FF (DSS) varying all three scales within = — 2p^ (the dashed lines indicate a 

±20% uncertainty). Middle: PDF uncertainty obtained for fixed ^ = p^ and FF (DSS) with three 
PDFs: CTEQ6.6, MSTW08, NNPDF1.2 (the dashed lines indicate ±8% differences). Bottom: FF 
uncertainty obtained for fixed scales {fi = p^) and PDF (CTEQ6.6) with three FFs: AKK08, DSS, 
HKNS (the dashed lines represent ±20%). _ -ik _ 



To assess the uncertainties linked to the choice of PDFs, FFs and scales // in the do- 
main of energies covered by the LHC, we have computed pp ^ X sd, a. fixed y/s = 7 TeV 
for various combinations of the theoretical ingredients as done for the Tevatron prediction 
(see Section Figure |^ (top) shows that the scale uncertainty is smaller (±20% above 
Pj, ~ 10 GeV/c) than found at lower (Tevatron) energies (see Fig. ^top). The middle plot 
of Fig. ^ shows that the uncertainty linked to the PDF choice is also slightly smaller than 
found at Tevatron, of the order of ±8%. Finally, the bottom panel shows that the fractional 
FF uncertainty is at most of ±20% above Pj, ~ 10 GeV/c, whereas below that transverse 
momentum the uncertainties increase up to ±40%. In the range Pj, ^ 30 - 60 GeV/c the 
FF choice has uncertainties of only 10 percent. Those results point again to a somehow 
smaller FF uncertainty than found at Tevatron (see Fig. ^bottom). A simple quadrature 
addition of the fractional uncertainties linked to the scales, PDF and FF choices results 
in a total theoretical uncertainty of around ±35% for the NLO single inclusive charged 
hadron spectrum in p-p collisions at LHC energies. 

4.2 Interpolation of measured charged- hadron spectra at ^/s = 5.5 TeV 

One of the assets of the successful RHIC physics program has been the ability to study the 
production of hard processes in p-p and nucleus-nucleus ( A-A) collisions at the same centre- 
of-mass energy. At the LHC, protons and ions have to travel in the same magnetic lattice^^ 
i.e. the two beams are required to have the same charge-to-mass ratio Z/A. This limits the 
beam momentum of a given species top = 7 TeV x Z/A for the nominal 8.3 T dipole bending 
field. The nominal nucleon-nucleon cm. energy for Pb-Pb collisions at the LHC is thus 
^ys^ = 5.5 TeV for lead ions with A = 208 and Z = 82. Since the maximum cm. energy 
in the first LHC p-p runs is half of the nominal value, ^/s = 7 TeV in lieu of y/s = 14 TeV, 
the first Pb-Pb runs are actually expected at a maximum ^/s^ = 2.76 TeV. In order to 
correctly normalize the yields measured in Pb-Pb collisions at y^s^ = 2.76,5.5 TeV, it 
will be crucial to get reliable estimates of the corresponding cross sections in p-p collisions 
at the same cm. energy. Ideally the predictions in p-p collisions should take advantage of 
the data accumulated at the LHC at nearby energies and be obtained with the smallest 
model-dependence possible to avoid any theoretical prejudice. In the following we present 
two methods for rescaling experimental p-p charged hadron spectra, measured at different 
cm. energies than those expected for heavy-ion collisions, based respectively on pQCD 
yield ratios and Xj, -scaling. 

4.2.1 Centre-of-mass energy rescaling 

As seen in the previous section, at a given cm. energy there are combined uncertainties of 
the order of ±35% in the NLO predictions for the absolute -differential cross sections of 
charged hadrons at LHC energies. Most of these uncertainties - in particular the largest 
scale dependence - however cancel out when taking ratios of the predicted perturbative 

^^The magnetic rigidity is defined as p/Z — Br for an ion witli momentum p and charge Z that would 
have a bending radius r in a magnetic field B. 
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yields at different, yet close, cm. energies. One can, thus, rescale the p-p spectrum 
measured at a given ^/s = X TeV (say, 7 TeV) to a = 2.76, 5.5 TeV reference value 
with a simple prescription: 

dcTref (^/i = 2.76, 5.5 TeV) ^ / da^^^/dp^ (^i = 2.76,5.5 TeV) \ ^ d^cxp = X TeV) 
dp^ V 'i-f^NLo / ^Pt iVs = X TeV) J dp^ 

(4.1) 

As an example, we plot in Fig. ^ the scaling factors as a function of Pj, obtained from the 
ratios of 2.76-TeV/7-TeV and 5.5-TeV/7-TeV pQCD yields. They have been obtained with 
up to 14 combinations of the different theoretical ingredients (scales, PDFs and FFs). The 
important scale dependence (Fig. 0, top) largely cancels out and only residual differences 
arise from the slightly different parton x and hadron z momentum fractions probed at 
different energies (of course, the closer the cm. energies the smaller the uncertainty in the 
yield ratios) . The maximum theoretical uncertainties of the rescaling factors amount to a 
small ±5% (resp. ±12%) for ^/s = 5.5 TeV (resp. 2.76 TeV), which will be likely below 
the expected uncertainties in the 7- TeV p-p experimental spectrum alone. 
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Figure 8: Rescaling factors of the -differential charged hadron cross-sections in p-p at 
-y/s 7 TeV down to lower y/s = 2.76, 5.5 TeV values, obtained from the ratio of the corresponding 
NLO calculations, Eq. ([4.l|) . The various curves show the small residual differences arising from 
different scale/PDF/FF choices. 
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4.2.2 Xy-scaling 



In this last section, we suggest to use the a:j,-scahng of particle production in high-energy 



scattering discussed in Sect. |3^ in order to predict the large-pj, hadron production cross 
sections in p-p collisions at ^/s = 2.76 and 5.5 TeV from the interpolation of Tevatron 
(Vi = 1.96 TeV) and LHC {^/s = 7 TeV) data. 

Assuming that Eq. ( f3.1| ) holds from Tevatron to LHC^^, it is straightforward to deduce 
the invariant cross section at a given and Xj, from previous measurements performed 
at Tevatron and LHC. Using the power-law interpolation^^, Eq. ( |3.1| ), the invariant cross 
section cr'"^ = E Sa/d?p reads 



cj""{l.m TeV,Xj 



^inv(7 TeV,Xy) 
o-i°^(1.96 TeV,Xj 



(4.2) 



where we define a = ln(y^/1.96)/ ln(7/1.96). The relative uncertainty on the cross section 
at ^/s resulting from the power-law interpolation from the data at 7 and 1.96 TeV is thus 
simply given by 



a" 



'{^/s■, 



'(1 



5cji'^^(1.96 TeV,x, 
(t'"^(1.96 TeV,Xj, 



5^inv(7 TeV,Xj,) 
^inv(7 TeV,Xj,) 



(4.3) 



where (5cr™^(1.96 TeV,Xy) and (5(T™^(7 TeV,Xj,) are the uncertainties of spectra mea- 
sured at Tevatron and LHC, respectively. Let us suppose for simplicity that the exper- 
imental relative uncertainty (5o""™/cj™^ on the hadron spectrum is identical at Tevatron 
and LHC, thus the relative uncertainty on the interpolated cross section at \fs will be 
a/o^ -h (1 - of X baja, i.e. 0.83 (^c7'"^/(t'"^ at both ^ = 2.76 and 5.5 TeV. Therefore, 
if the measurements at Tevatron and LHC used for the interpolation are precise enough, 
this procedure would allow for predictions whose uncertainties become possibly smaller 
than the usual theoretical uncertainties of NLO QCD calculations. Note however that the 
uncertainty Eq. ( fl.3D only reflects the propagation of errors in the power-law interpolation 
and does not account for the systematic uncertainty of the procedure currently used. 

Since the invariant cross sections are compared at a given x^,, the prj, range reached 
at the lower (Tevatron) and upper (LHC) limits of the interpolation domain is crucial. 
The currently "reliable" Tevatron data extend up to ~ 20 GeV/c which allows for a 
prediction up to Pj, ~ 30 (60) GeV/c at -^/i = 2.76 (5.5) TeV. Conversely, the upper limit 
at i/i = 2.76 (5.5) TeV is 40% (80%) of the highest p^ to be reached at i/i = 7 TeV. 

In order to check this procedure, the Pj,-spectrum of mid-rapidity charged hadrons at 
\fs = 5.5 TeV has been estimated from the pythia spectra at ^/s = 1.96 and at 7 TeV 
using the interpolation Eq. ( [4. 2] ) with a = 0.81, within the range Xj, w 2-10-3 - 0.2. The 



^^We note that the Tevatron data are measured in p-p coUisions unhke the p-p coUisions at the LHC. The 
differences between both systems on unidentified hadron production at midrapidity is very smali (especialiy, 
far away from the valence quark region, i.e. for x^, ^1). 

^^Note that the exponent n depends in principle on p^ (and thus Xj,) from the scaling violations in QCD 
and should approach n = 4 in the Bjorken limit. However this dependence is expected to be logarithmic 
and can be safely neglected in the p^-range being considered here. 
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result is plotted in Fig. |^ and compared to the direct pythia calculation of the hadron 
spectrum at y/s = 5.5 TeV. As it can be seen, the x^- interpolated cross section reproduces 
nicely the MC result at = 5.5 TeV above « 5 GeV/c. 




5. Summary 

We have compared the latest high-p^, charged particle spectrum measured by CDF in 
proton-antiproton collisions at ^/s = 1.96 TeV to various perturbative QCD expecta- 
tions based on next-to-leading-order calculations (incnlo), parton-shower Monte Carlo 
(pythia), and Xy-scaling respectively. The NLO calculations employ the latest sets of 
parton distribution functions (PDFs) and fragmentation functions (FFs). The Tevatron 
data can be well reproduced below p^. 20 GeV/c for the choice of scales ^ = 2pj., 
CTEQ6.6 parton densities, and AKK08 parton-to-hadron fragmentation functions. Above 
this value, the CDF spectrum starts to rapidly deviate, by up to three orders of mag- 
nitude, from the predictions. The most important source of theoretical uncertainty, of 
order ±30%, is related to the choice of the factorization, fragmentation and normalization 
scales. The maximum uncertainties linked to the choice of the PDFs and FFs are ±10% 
and ±25% respectively. A conservative quadratic addition of all these differences results in 
a maximum ±40% uncertainty in the NLO calculations which cannot by any means explain 
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the important data-theory disagreement above k, 20 GeV/c. 

We have next determined with the PYTHIA MC the possible extra contributions of 
high-pj, charged particles, including leptons, coming from heavy-quark fragmentation as 
well as from real and virtual vector-boson production either single-inclusive or in associa- 
tion with a jet. The addition of such processes, which amount to about an additional ten 
percent of the charged particle yield above ~ 40 GeV/c, does not help to reduce the 
large data-theory deviation. The CDF spectrum also fails to fulfill simple scaling expec- 
tations which are empirically confirmed by all other high-p^, hadron spectra measured so 
far in p-p collisions in the range -v/i = 0.2 - 1.8 TeV. Moreover, the power- law exponent of 
the CDF data above x^, k, 0.02, is below the n = 4 limit expected from simple dimensional 
arguments for pure 2^2 parton scattering in QCD. 

We conclude that the fact that the NLO calculations largely fail to reproduce the mea- 
sured CDF single hadron spectrum at large p^, while simultaneously reproducing correctly 
the single jet -differential cross sections, and that the measurement violates simple phe- 
nomenological expectations such as x^-scaling, point to a possible experimental problem 
in the Tevatron data above Pj, ~ 20 GeV/c (or to unknown sources of charged particles 
not considered here, a possibility disfavoured in [43]). 

The NLO predictions of charged hadron spectra at LHC energies -y/s = 0.9-14 TeV, 
have also been provided. Finally, we have proposed two simple interpolation procedures, 
based on a pQCD-rescaling and an x^-scaling of (future) experimental p-p transverse mo- 
mentum spectra, in order to obtain the nuclear modification factors of high-p^, charged 
hadron production in nucleus-nucleus collisions at intermediate {^/s]^ = 2.76, 5.5 TeV) 
LHC energies. 
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